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ABSTRACT 

As a first step for studying star formation in tlic extreme outer Galaxy (EOG), we obtained deep 
near-infrared (J, H, if -bands) images of two embedded clusters at the northern and southern CO 
peaks of Cloud 2, which is one of the most distant star forming regions in the outer Galaxy (galactic 
radius Rg - 19 kpc). With high spatial resolution (FWHM - 0'.'3~0'.'35) and deep imaging {K - 21 
mag, 5 cr) with the IRCS imager at the Subaru telescope, we detected cluster members with a mass 
detection limit of < O.IM©, which is well into the substellar regime. These high quality data enables 
a comparison of EOG to those in the solar neighborhood on the same basis for the first time. Before 
interpreting the photometric result, we have first constructed the NIR color-color diagram (dwarf 
star track, classical T Tauri star (CTTS) locus, reddening law) in the Mauna Kea Observatory filter 
system and also for the low metallicity environment since the metallicity in EOG is much lower than 
those in the solar neighborhood. The estimated stellar density suggests that an "isolated type" star 
formation is ongoing in Cloud 2-N, while a "cluster type" star formation is ongoing in Cloud 2-S. 
Despite the difference of the star formation mode, other characteristics of the two clusters are found 
to be almost identical: (1) K-haiid luminosity function (KLF) of the two clusters are quite similar, as 
is the estimated IMF and ages (~ 0.5-1 Myr) from the KLF fitting, (2) the estimated star formation 
efficiencies (SFEs) for both clusters are typical compared to those of embedded clusters in the solar 
neighborhood (~ 10 %). The similarity of two independent clusters with a large separation 25 pc) 
strongly suggest that their star formation activities were triggered by the same mechanism, probably 
the supernova remnant (GSH 138-01-94) as suggested in Kobayashi et al. and Yasui et al. 
Subject headings: infrared: stars ~ stars: formation ~ stars: pre-main-sequence - open clusters and 

associations: individual (Digel Cloud 2) - stars: luminosity function, mass function 

- ISM: clouds - supernova remnants 



1. INTRODUCTION 

The extreme outer Galaxy (hereafter "EOG"), de- 
fined here as the region at Galactic radius (Rg) greater 
than 18 kpc, could be regarded as the outermost re- 
gion of the Galaxy because the distribution limits 
of Population I and II stars are known to be 18 - 
19 kpc and 14 kpc, respectively (jPigel et all I1994D . 
The EOG is known to h a ve very low gas d ensity 
(|Brand fc WouterlootI [19951: iNakanishi & Sofud \mm 
and low metallicitv (jSmartt fc Rollestonlll997D . Such an 
extreme environment could serve as a "laboratory" for 
studying the star formation process in a totally different 
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environment from that in the solar neighborhood. The 
star-forming process in such environments has not been 
studied extensively as for the nearby star-forming regions 
because of the large distance. The EOG does not seem 
to have long and complex star formation history because 
there is little or no perturbation from the spiral arms, 
which is the major and powerful trigger for continuously 
converting a large amount of molecular gas into stars. 
Therefore, in the EOG there is a good chance to observe 
the detail of the each star formation process, such as 
the triggering process of cluster formation, without be- 
ing disturbed by the entangled star formation history in 
space and time. 

There are some synthetic studies of molecular clouds 
in EOG and far outer Galaxy (hereafter "FOG"), whose 
Rg is not more than 18 kpc but no less than 15 kpc. 
The pioneering survey of molecular clouds and star 
formation beyond the solar circle was undertaken by 
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IWouterloot fc Brandl (|1989f) . who identified the ^^CO 
molecular emission associated with color-selected IRAS 
point sources and used kinematic distances to place 
the sources within the Galaxy. Many of the IRAS 
sources were found to lie in the far outer Galaxy. 
Follow-up studies of the gas properties of many of 
th ese far outer Galaxy molecul a r clou ds were presented 
in iBrand fc Wouterlooi] ([1994 I1995D se parately. In 
an independent search. iDigel et al.l ( 1994[ ) found eleven 
molecular clouds in the EOG, based on CO obser- 
vations of di stant H I peaks in the Maryla nd- Green 
Bank survey l|Westerhout fc Wendlanp 11983 ). Follow- 



hensivc study of a star forming region, WB89-789, whose 
distance is estimated at D = 11.9 kpc, Rn ^ 20.2 kpc. 
This region was found in lWouterloot fc Brandl (|1989[ ) de- 
scribed above. They obtained NIR imaging with limiting 
magnitude of mx = 16.5 and 17.5 mag for 10 a and 5 a, 
respectively, which correspond to a 2-3 Mq dwarf main 
sequence star with 10 mag of visual extinction. They 
also used molecular line and dust continuum observa- 
tions for properties of the molecular cloud and showing 
the existence of an outflow. They also obtained optical 
spectroscopy o f a star for c o nfirm ation of the distance 

, to the region. ISantos et al.l (|2000( l discovered two dis- 

ing those pioneering works, iHever et al.l (|1998l ) con- tant embedded young clusters in the FOG {Rg = 16.5 

kpc; D = 10.2 kpc) with near-infrared imaging of CO 
clouds associated with IRAS sources. Their detection 
limit {iriK — 16.4 mag) corresponds to a 1 Myr old star 
of ~l-2 Mp) seen through 10 mag of visual extinction. 
ISnell et all 1)2002! ) have conducted a comprehensive NIR 
survey of embedded clusters in the FOG ( Rp = 13.5- 
17.3 kpc) based on the FCRAO CO survey (jHever et all 
Il998f ). They found 11 embedded clusters with the de- 
tection limit of mii-'=17.5 mag. For their most distant 
cluster, this magnitude corresponds to a 1 Myr old, 0.6 
Mq star with no extinction, or a 1 Myr old, 1.3 Mq star 
seen through 10 mag of visual extinction. These indepen- 
dent ncar-infrarcd searches so far found cluster members 
with the mass detection limit of ~ 1 Mq. 

As a next step to study the details of the star forma- 
tion activity in the EOG, we obtained deep and high 
spatial resolution NIR images of the two most distant 
embedded clusters in Digel Cloud 2 with Subaru 8.2-m 
telescope. Because of the achieved sensitivity to detect 
substellar mass members, we could compare embedded 
clusters in EOG to those in the sol ar neighborhood for 
the first time. In I Yasui et al.l ()2006l . hereafter Paper I), 
we presented the data for the Cloud 2-N cluster and (i) 
estimated the age of the cluster at less than 1 Myr by 
comparing the observed KLF to model KLFs based on 
typical IMFs, and (ii) suggested that the formation of the 
cluster was triggered by the SNR shell based on the age 
and the geometry of the cluster. In this paper, we present 
data of the Cloud 2-S cluster and also re-analyze the data 
of Cloud 2-N cluster to investigate the fainter stellar pop- 
ulation. For interpreting the photometric result, we have 
constructed for the first time the NIR color-color diagram 
(dwarf track, CTTS locus, and reddening vector) in the 
Mau na Kea Observatory (MKO ) near-infrared filter sys - 
tem (jSimons fc Tokunagal l2002l : iTokunaga et~an l2002f l. 
To accommodate to the low-metallicity in EOG, we also 
considered the mass-luminosity relation and the color- 
color diagram for low metallicity stars. The established 
methods here would be useful for future study of star for- 
mation in EOG and FOG, and ultimately nearby dwarf 
galaxies. We then discuss the KLF and star formation 
efficiency in the EOG based on the results of these two 
clusters. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Subaru IRCS imaging 

We obtained J (1.25 //m), H (1.65 //m), K (2.2 //m)- 
band deep images of the Cloud 2-N and -S clusters with 
the Subaru Infrared Camera and Spectrograph (IRCS; 
Tokunaga et al. 1998; Kobayashi et al. 2000; Ter- 
ada et al. 2004) with a pixel scale of 0'.'058 pixel"^. 



ducted a comprehensive CO survey of molecular clouds 
in the FOG with Five College Radio Astronomy Obser- 
vatory (FC RAO) CO survey. W it h Hq and radio contin- 
uum studv. iFich fc Blitd 11981 ■ iBrand fc Bliti (Il99l . 
and iRudolph et al.l ( 19961 ) identified a number of H II 
regions associated with the FOG clouds, indicating the 
presence of recent massive star formation. 

For the study of star formation, it is ultimately nec- 
essary to observe YSOs (young stellar objects) at near 
infrar ed (NIR) wavelengths. iKobavashi fc Tokunagal 
()2000f ) found YSOs associated with a molecular cloud 
in EOG for the first time with a wide-field NIR survey. 
This c loud is denoted as Cloud 2 in the list of lDigel et al.l 
(fl99l and is one of the most distant molecular clouds 
in the outer Galaxy region. It is located at a very large 
Galactic radius {Rg = 18-28 kpc) in the second quad- 
rant of the Galaxy (/ ~ 137.75°, b 1.0). It is by far 

the largest molecular cloud among those found in EOG 
with a mass comparable to that giant molecular clouds 
(M ~ 2 X lO'' Mq] Digel et al. 1994). The Galactic 
radius of Cloud 2 has been estimated by various meth- 
ods at Rg = 28 kpc (heliocentric distance: D = 21 kpc) 
from the kinematic distance of the CO cloud (jPigel et al.l 
119941 ): Rg = 23.6 kpc (D = 16.6 kpc) from the latest H I 
observation (jStil fc IrwinI l2001l) : and Rg = 15-19 kpc 
[D = 8.2-12 kpc) from the detailed o ptica l study of a B- 
type star, MRl (jMuzzio fc Rvdgrenlll973 : 'Smart t et all 



19961), which is associated with Cloud 2 (dc Gcu s et al.l 
19931 ). Throughout this paper we adopt i?g = 19 kpc 
(D = 12 kpc) because it is about the mean value of all the 
estimated distances and because spectroscopic distance 
of stars should be more accurate than kinematic distance. 
Metallicity at the Galactic radius of Cloud 2 is estimated 
to be ^ — 1 dex from the standard metallicity curve by 
ISmartt fc RoUestonl (|1997l l. In fact, the metallicity of the 
B-type star, MRl, is measured at — 0.5 to —0.8 dex by 
optical high-resolutio n spectroscopy (jSmartt et al.l[T996l : 
iRolleston eFall l2000l ) and that of the Cloud 2 itself is 
meas ured at —0.7 dex from the radio molecular emission 
lines ([Lubowich et al.l l2004l : ?). These metallicity val- 
ues are comparable to that of LMC (~ —0.5 dex) and 

SMC ( 0.9 dex: lArnault et al.lll988l ). Using deeper 

near-infrared imaging with QUIRC at the University of 
Hawaii (UH) 2.2m telescope, Kobayashi et al. (2007) 
also found two embedded clusters at the northern and 
southern CO peaks of Cloud 2. They suggested that the 
overall star formation activity in Cloud 2 was triggered 
by the supernova remnant (SNR) H I shell, GSH 138-01- 
94, and that supernova triggered star formation could be 
one of the major star formation modes in EOG. 
Recently Brand fc Wouterloot (2007) reported compre- 
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IRCS employs the Mauna Ke a Observatory (MKO) near- 
infrared photometric filters (jSimons fc Tokunagal 120021 : 
iTokunaga et al.ll2002[ ). The entire infrared cluster was 
sufficiently covered with the ~ I' field of view. For 
the Cloud 2-N cluster, the observation was conducted 
on 2000 December 1 UT with the total integration time 
of 600, 450, 900 s for J, H, and K bands, respectively. 
For the Cloud 2-S cluster, the observation was conducted 
on 2000 October 25 UT with the total integration time of 
600, 600, 420 s for J, iJ, and K bands, respectively. The 
observing condition was photometric and the seeing was 
excellent (~ 0'.'3-0'.'35) throughout the observing period 
on both nights. 

2.2. Data Reduction 

All the data for each band were reduced with IRAF^ 
with standard procedures: dark subtraction, flat field- 
ing, bad-pixel correction, median-sky subtraction, image 
shifts with d ithering offsets, and combining. In our pre- 
vious paper ijPaper II ) , we used the DAOFIND algorithm 
in IRAF DAOPHOT package with the detection thresh- 
old of 10 cr for the source detection in the Cloud 2-N 
cluster. In the present study, we performed source detec- 
tion for the Cloud 2-N cluster again and the Cloud 2-S 
cluster using SExtractor (Bertin et al. 1996) with the 
detection threshold of 5 cr above the local background. 
Photometry was performed using IRAF APPHOT pack- 
age with aperture diameters of 0'.'58 (10 pixel) and 0'.'46 
(8 pixel), for Cloud 2-N and -S, respectively. The aper- 
ture sizes were chosen to achieve the highest S/N. For 
four pairs of binary stars in the Cloud 2-N cluster, 0'.'35 
(6 pix) diameter aperture was applied with an aperture 
correction because the separation is less than 0'.'5. For 
photometric calibration we used four and seven bright 
and isolated stars for the Cloud 2-N and -S cluster, re- 
spectively, in the images for which photometry was ac- 
curately done later with newly obtained JHK images 
of the fields with Subaru's new wide-field near-infrared 
camera MOIRCS (Ichikawa et al. 2006). The flux un- 
certainties are estimated as the same way as in Paper I. 
The resultant limiting magnitudes (5 cr) for the Cloud 
2-N cluster are J = 22.0, H = 21.0, and K = 20.9, while 
those for the Cloud 2-S cluster are J = 21.8, H = 21.0, 
and K = 20.4. The pseudocolor images of the observed 
field and the spatial distributions of detected stars arc 
shown in Figs. 1 and 2 for the Cloud 2-N cluster and -S 
cluster, respectively. 

3. COLOR-COLOR DIAGRAM IN MKO PHOTOMETRIC 
SYSTEM 

For identifying of young cluster members, JHK 
color-color diagram with dwarf and giant star tracks 
in Johnson-Glass system (iBessel &: Brettj 1 19881 here- 
after "BfcB"), CTTS locus in CIT system (IMever et al.l 
Il997t). and r eddening vector in Arizona system 
(iRieke fc Lebofs kv 198^ have been widely used (e.g., 
iLada et all 12000) . Because the Mauna Kea Observa- 
tory (MKO) filter sy stem (jSimons fc Tokunagal [200^ : 
iTokunaga et al.l r2002f ). which we used for the present 
study, is relatively new, the color-color diagram in this 

^ IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which arc operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



system have not yet been clearly established. Although 
we used the color-color diagram with t he dwarf track etc. 
in other systems in our previous paper (jPaper J ), we con- 
structed a color-color diagram solely in the MKO system 
for more rigorous treatment of the data. 

Firstly, we plotted colors of dwarf stars^ observed with 
MKO filters (Hewett et al. 2006, open circle, B9V- 
M6V^ Leggett et al. 2002, plus, M4V-M6V''; Leggett 
et al. 2006, cross, M4-M6^) to construct a dwarf track 
in the MKO system (solid line). Although this MKO 
track appears to be almost identical to the B&B track 
for stars earlier than K type stars, the MKO track has 
a smaller J — H color than the B&B track for M type 
stars. Because the MKO J filter has a narrower band 
width compared to the ot her filter systems (see Fig. 5 in 
[Stephens fc Leggett|[200l . the MKO J magnitude is less 
affected by the strong II2O vapor absorption bands of M 
type stars. This should result in a brighter J magnitude, 
thus the smaller J — H color in the MKO system. We 
interpret this as the reason for the shift of the track po- 
sition for M type stars and also the reason why the shift 
becomes larger with later spectral type. 

Secondly, in order to construct the CTTS locus in 
the MKO system, we checked the 2MASS J, H, /v-band 
magnitudes of classical T Tauri stars® in Strom et al. 
(1989), which were use d to derive the orig inal CTTS lo- 



cus in the CIT system iMever et al.l ()1997h . After dcred- 
dening th e 2MASS co l ors of these stars following the 
method of IMever et al.l (|1997f ) with the reddening law in 
iCambresv et al.l ( 20021 ) . we converted the 2MAS S col ors 
to the MKO colors using the conversion in .Leggett et al] 
()2006f) . We derived the CTTS locus in MKO system 
by least-square fitting of the MKO colors: the resultant 
locus is (J - iJ)MKO = (0.32 ± 0.07)(i7 - K)mko + 
(0.54 ± 0.04). Compared to the original CTTS locus. 
{J-H)ciT = (0.58±0.11) X (i7-is:)ciT + (0.52 ±0.06), 
the MKO locus has smaller J — H with increasing H — K. 
This tendency is consistent with that for the dwarf track 
as discussed above. 

Lastly, we considered the reddening vector, which is 
the the slope (J — H)/{H — K) and the length of the 
Ay = 5 mag vector. Using the near -infrared camera 
SIRIU S equipped with the MKO filters. iNishivama et al.l 
(|2006f) estimated that {J-H)/{H-Ks) = 1.72 from the 
data toward the Galactic Center. This value is almost 
identical to the con ventional slope (J — H /H — K = 1.70) 
in Arizona system (jRieke fc Lebofskv! Il985l ). Although 
IRSF uses a Ks band filter instead of a K band filter, 
th e difference of the filt er profile is so small (see Fig. 1 
in ITokunaga et aLll2002f) that the magnitude differences 
are negligible^. Therefore, we could conclude that the 

^ Plotted are only for dwarfs earlier than M6, which is the end 
point of the B&B dwarf track. 

^ I Hewett "et al., (.2008 ) shows B-Y dwarf stars. These are com- 
puted from their own spectra, (synthetic colors) 

* ILeggett etlj] 12003) shows the colors of M4V-T8V stars. 

^ ILeggett et al.l H2006I ') shows the colors of B7.5-M9.5 stars. We 
plotted dwarf stars and also stars without clear designation of lu- 
minosity class. 

^ Although IMever et all l(T997h used 30 CTTSs from Strom ct 
al. (1989) for their derivation, we found 33 usable CTTSs with the 
spectral type range of G0-M5 as defined in Strom et al. (1989). 
Wc removed four stars from our analysis because their 2MASS 
photometric quality is not good. 

^ For reference, we compared K and Ks magnitudes of Persson 
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reddening veetor is almost identical for both Arizona and 
MKO systems. In the following section, we use the MKO 
color-color diagram for all the analysis. 

4. COLOR-COLOR DIAGRAM FOR LOW METALLICITY 
ENVIRONMENT 

Because the metallicity in EOG is significantly lower 
than in the solar neighborhood, we should consider the 
effect of low metallicity o n the color-color diagram. Fol- 
lowin g the PMS models bv lD'Antona fc Mazzitellil () 19971 . 
fT998l) . we suggested that NIR absolute magnitudes of 
PMS stars with —0.3 dex (1/2 solar metallicity) is not 
sign ificantly different from those with the solar metallic- 
ity (jPaper J ). Therefore colors also should not signifi- 
cantly change down to a metallicity of —0.3 dex. In this 
section we discuss the effect of low metallicity extended 
to —1 dex on the J — H and H — K colors considering 
b oth obse r vation s and models in the literature. 

iLeggettl (|1992l ) presented NIR photometry of low-mass 
stars in young disk ([M/H] ^ 0.0), young-old disk, old 

disk ([M/H] 0.5), old disk-halo ([M/H] 1.0), 

and halo ([M/H] < —1.0), showing that J — H colors 
of dwarfs between K7 (0.5 Mq) and M5 (~ 0.1 Mq) 
types decrease by ^ 0.1 mag with decreasing metallic- 
ity (Fig, m. Simila r tendency is seen in the model by 
iGirardi et all (|2002t ). who listed the NIR magnitudes of 
stars with mass of 0.15-1 Mq for various metallicities 
(Z = 0.001,0.004,0.008, and 0.019). Although the color 
values are slightly different for the observation and the 
Girardi model, both show the quite similar color varia- 
tion A{ J-H) = -0.1 for A [M/H] 1, which strength- 
ens the validity of the metallicity dependence. However, 

for the metallicity of Cloud 2 ([M/H] 0.5), the dwarf 

track is still almost identical to the standard B & B dwarf 
track within the uncert ainty (Fig. [H). As fo r the stars 
with mass of > 1 Mq, IGirardi et al.l (|2002l ) suggested 
t hat the colors wea.kly d epend on metallicity. 

iNakaiima et all ()2005[ ). who observed the LMC (-0.5 
dex) in the MKO system, suggested that the slope of ( J— 
H)/{H — K) for LMC appears to be consistent with that 
for solar metallicity (Nishiyama ct al. 200_1); suggesting 
that the reddening vector for the LMC metallicity is not 
significantly different from that for the solar metallicity. 

We conclude that dwarf track and reddening vector 
in JHK color-color diagrams do not strongly depend 
on the metallicity down to that of the stars in Cloud 
2 (^ —0.7 dex). Therefore, we do not include the ef- 
fects of low metallicity in the following analysis using 
color-color diagram, such as the identification of cluster 
members (§ 15. ip . the estimation of extinction and color 
excess (§ I5.2|1 . 

5. RESULTS 

5.1. Identification of Cluster Members 

We conducted identification of the cluster members us- 
ing pseudocolor pictures of the observed fields (Fig. [TJ [2]) 
and (J — H) vs. {H — K) color-color diagrams of all the 

standard stars IjPersson et al.|[l998f ) and found that the magnitude 
differences are almost within 0.02 mag (the average difference is 
— 0.0070±0.0067 ma g). Although these m agnitudes are in the LCO 
photometric system IjPersson et ahiriggSI I. similar results would be 
expected for the MKO photometric system because the difference 
of the MKO Ks and K filter profiles is smaller than that of LCO 
Ks and K filter profiles. 



detected sources (Fig. [5]). Wc identified cluster members 
which have color oi {H — K) > 0.5 and are within the 
cluster regions as identified in Kobayashi et al. (2007). 
Because there are no foreground molecular clouds in front 
of Cloud 2 (Kobayashi & Tokunaga 2000), the red colors 
of the cluster members should originate from the extinc- 
tion by Cloud 2 itself as well as by circumstellar material 
(see § 5.2 for more discussion). Contamination from the 
background sources should be very small in view of the 
large Rg of Cloud 2. 

We found 72 cluster members in the Cloud 2-N clus- 
ter and 66 cluster members in the Cloud 2-S cluster. 
We also identified 10 YSOs in the Cloud 2-S sub-cluster 
(Kobayashi et al. 2007), which is ~ 40" northern-east 
from the Cloud 2-S cluster. Fig. [1] and Fig. [2] show the 
spatial distributions of these cluster members and field 
stars in the Cloud 2-N and -S fields, respectively. Be- 
cause of the re-analysis with the fainter detection limit, 
the number of identi fied Clou d 2-N cluster members has 
increased from 52 in iPaperl to 72. 

5.2. Extinction and Disk Color Excess of Cluster 
Members 

The extinction and disk color excess for each star were 
estimated using the color-color diagram. For reliable es- 
timation of the parameters, only stars in the positions 
reddened from the CTTS locus and the dwarf locus are 
used (42 cluster members and 36 field stars for the Cloud 
2-N cluster, 45 cluster members and 45 field stars for the 
Cloud 2-S cluster). For convenience the dwarf locus was 
approximated by the extension of the CTTS locus drawn 
out to H — K ~ 0.1 mag. In the color-color diagram 
the extinction Ay of each star was estimated from the 
distance along the reddening vector between its location 
and the stellar loci. The resultant Ay distributions of 
Cloud 2-N and Cloud 2-S cluster members (Fig. [6]) have 
peak at 7.2 and 6.1 mag. respectively {Ak = 0.81 and 
0.68 mag) and those for field stars in the both Cloud 2 
fields have a peak at 2.2 mag {Ak = 0.25 mag). The dis- 
tributions of the both clusters' members show that the 
both clusters are reddened uniformly by about Ay ~ 4-5 
mag inside Cloud 2. This verifies the selection method 
of cloud members as described in § 15.11 These values of 
Ay are also consistent with the estimated from sub radio 
and continuum observation (Ruffle et al. 2007). 

We also constructed the distributions of the unred- 
dened color {H—K)o of the cluster members and the field 
stars (Fig. [7]). For both clusters, the distributions of the 
cloud members and field stars show the clear peak offset 
of about 0.1 mag. The difference of the average {H — K)o 
between the cluster members and the field stars (Cloud 
2-N cluster: 0.26 mag. Cloud 2-S cluster: 0.37 mag) can 
be attributed to thermal emission of circumstellar disks 
of the cluster members. Assuming that disk emission 
appears in the K band but not in the H band, the disk 
color excess of the Cloud 2-N and -S cluster members 
in the K band, AiCdisk, is equal to 0.26 and 0.37 mag, 
respectively. 

6. DISCUSSION 

6.1. Stellar Density Variation 

iLada fc Ladal ()2003l ) compiled a catalog of embedded 
clusters that are located in 13 < 2 kpc (distance modu- 
lus DM < 11.5 mag) and detected cloud members with 
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the apparent lim iting magnitudes of tuk ~ 16.0 mag. 
ICarpented ()2000f ) analyzed 2MASS data of the Perseus, 
Orion A, Orion B, and Monoceros R2 clusters at I? ~ 500 
pc (DM 2± 8.5 mag) and detected members with apparent 
limiting magnitudes of nixs ~ 14.3 mag. Both catalogs 
thus consist of clusters that are detected with absolute 
limiting magnitude of Mk ~ 5-6 mag in the solar neigh- 
borhood. In the present study, since the Cloud 2 clus- 
ters at D ~ 12 kpc (DM = 15.4 mag) are detected with 
apparent magnitudes of niK ^ 20.5 mag, the absolute 
limiting magnitudes is Mk 5 mag. Therefore, we can 
compare the stellar density of embedded clusters in EOG 
to those of embedded clusters in the solar neighborhood 
o n the same ba s is for the first time. 

I Adams et al.l ()2006l) found a correlation between the 
number of stars in a cluster and the radius of the cluster, 
using the tabulated cluster properti es in iLada fc Ladal 
((200a) and lCaroente^ (pOOl (Fig.^. lAllen et al.l ((2007f ) 
suggested that the average stellar density of clusters have 
a few to a thousand members varies by a factor of only 
a few. Stellar densities of the Cloud 2-N and -S clusters 
were estimated at 13 pc^^ and 50 pc^^ from the spa- 
tial distribution of the identified cluster members in the 
2x2 pc^ area and a circle of 0.6 pc radius, respectively. 
For the Cloud 2-N cluster, the a bove den sity is a lit- 
tle larger than our old estimate in iPaperl (~ 10 pc~^) 
because we included the newly identified fainter stars. 
In Fig. [51 the Cloud 2-N cluster appears to be a kind 
of "loose cluster", while the Cloud 2-S cluster appears 
to be a kind of "dense cluster". For the typical stellar 
density of the clusters (~ 50 pc~^, see Figure the 
average projected distance between two adjacent stars is 
~ 0.14 pc, which corresponds to 2" at Z? = 12 kpc. 
For even densest clusters with 1000 pc^^, the average 
projected distance is ~ 0.03 pc, or 0'.'5 at = 12 pc. 
These separations are sufficiently resolved for this obser- 
vation. Even if all the stars are unresolved binaries, the 
"loose" and "dense" characteristics of the Cloud 2 clus- 
ters do not change (see gray circles in Figure [8]). The 
difference of the stellar densities is very interesting, be- 
cause the two clusters show very similar characteristics 
for other points as shown in the following subsections. 
Kobayashi et al. (2007) attributed the difference of the 
star formation mode to the ambient pressure difference 
of the two clusters. 

6.2. IMF and Age of the Clusters 

if-band luminosity functions (KLFs), which are the 
number of stars as a function of K-hand magnitude, of 
different ages are known to have different peak mag- 
nitudes and slopes (Muench et al. 2000). IMFs and 
ages of young clusters can be estimated by comparison 
of the observed KLF with model KLFs, which are con- 
structed from IMFs and M-L relat ions. I n this sec- 
tion, we use the model KLF from iPaper II (see § 4). 
Although the popular PMS evolut i on models s uch a s 
iD'Antona fc Mazzitellil(fT997l . fTool . iSiess et al.l 12000), 
and iBaraffe et al.l (llQQTT l. do not have models for metal- 



hcity less than 1/2 solar (—0.3 dex). We assumed that 
the M-L relation for solar metallicity can be applied to 
model KLFs for low metallicity because the difference of 
M-L relations for the solar metallicity and those for 1 /2 
solar metallicity is very little (see more detail in Paper 
I)- 



Following the method in lPaper^ . we constructed KLFs 
of the Cloud 2 cluster members (Fig. [9]). We estimated 
the detection completeness in each magnitude bin by 
putting artificial stars on random positions in the field 
and checking whether they are detected in the same way 
as for the real objects. Five stars are placed at a time 
in each magnitude bin and the check was conducted 200 
times, resulting in about 1000 artificial stars per magni- 
tude bin. The result is also shown in Fig. [9] (bottom), 
and the completeness-corrected KLFs are shown in Fig.[S] 
(top) with thick lines. Both KLFs appear to be very sim- 
ilar in the slope at brighter magnitude side and the peak 
position (~ 19 mag). 

The comparison of the completeness corrected KLFs 
with the model KLFs are shown in Fig. [TOl Because 
the slope of model KLF does not significantly depend 
on small d ifference of underlying IMFs (see Fig. 6, 7 
in iPaper J ), we used the Trapezium IMF ()Muench et all 
I2002t ). which is thought to be the most reliable IMF for 
young clusters (e.g., Lada & Lada 2003). A quick visual 
check of Fig. [TU] suggests that the ages of the Cloud 2-N 
and -S clusters are similar and it is ~ 0.5-1 Myr. It is 
difficult to estimate the age of the cluster with an accu- 
racy of 0.1 Myr because isochrone models f or ages of less 
than 1 Myr is thought to be uncertain fe.g.. IBaraffe et aTl 
I2002t l. However, we can at least conclude that the ages 
of the tw o clusters are no more than 1 Myr (see § 5.1 in 
IPaper t ). For the Cloud 2-N cluster, twenty cluster mem- 
bers are detected in the K ^ 20 magnitude bin and more 
in the fainter magnitude bins, while only four members 
were detected in the ~ 20 magnitude in (the faintest 
bin) iPaper J . However, this change in the fainter bins 
{K > 20) did not affect the age estimation, because the 
brighter bins (A' < 19) are more sensitive for the KLF 
fitting. The relatively good fit of the model KLF to the 
observed KLF suggests that the IMF at M > 0.1 Mg in 
Cloud 2 is not significantly different from the "universal" 
IMF for Trapezium (see also discussion in Paper I). 

The similarities of the KLFs and the estimated ages 
for the two independent clusters with a large separation 
(~ 25 pc) strongly suggest that their star formation ac- 
tivities were triggered by the same mechanism. The ages 
of the clusters were estimated at ~ 0.5-1 Myr, which is 
much less than that of the SNR H I shell (GSH 138- 
01-94; Stil fc Irwin 2001), 4.3 Myr. An independent 
age estimate of the Cloud 2-N cluster (~ 0.5-1 Myr) 
by Kobayashi et al. (2007) is in quite good agreement 
with our estimate. These age estimates strongly support 
the idea that the star formation activity in Clo ud 2 was 
triggered by the SNR H I shell, as discussed in IPaper II 
and Kobayashi et al. (2007). 

6.3. Star Formation Efficiency (SFE) 

The star formation efficiency (SFE), which is generally 
defined as [SFE = Mstars/(-/W"gas + M^tars)] where Mgas is 
cluster-forming core mass and Mstars is the total stel- 
lar mass, is one of the most fundamental parameters of 
the s tar and cluster-formation processes (jLada fc Ladal 
|2003[ ). The estimated SFEs of the Cloud 2 clusters could 
be "pure" SFE of one-time cluster formation event be- 
cause there has been no complex star formation history 
in the EOG. As discussed in § 16.11 the detection of the 
embedded cluster members in the EOG with the high- 
sensitivity {rriK — 21.0 mag) enables a fair comparison 
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of SFE in the EOG with that in the solar neighbor hood. 
However, because SFEs have been estim ated (e.g.. iLadal 
fl99l iLada et al.lfT997l : lHartmannl[200l) m various ways, 
ex. Mgas is estimated using different kinds of molecular 
emission lines (^^CO, ^^CO, and C^^O), Afstars is esti- 
mated with a very simple assumption that all stars have 
a mass of 1 Mq, 0.5 Mq, etc., here we attempted to es- 
timate SFEs of various star forming regions in a uniform 
fashion as described in the following paragraphs. 

For the estimates of Mgas we use ^'^CO because of 
the relatively high sensitivity even for distant molecu- 
lar clouds and also the wide availability of the archival 
data. Because Mgas from ^'^CO correlates well with those 
from C^^O (C^^O-i^CO mass correlation; iRidge et all 



120031 Fig. 18 (a)), it is possible to use C^^O data when 
^■^CO data is unavailable. For the accurate estimation 
of Afstars, it is ncccssary to take the mass of the low 
mass stars into account. We use a summary of Mstars 
for star forming mo lecular clouds in th e solar neighbor- 
hood {D <2 kpc: ILada fcTadal[200l Tab. 1). They 
derived Mstars for each cluster by predicting infrared 
source counts, down to 0.017 Mq stars, as a function 
of differing limit i ng ma gnitudes using the KLF models 
of iMuench et al.l ()2002D . As an example, we estimated 
the SFE of the well-known embedded cluster, NGC 1333, 
which is in the star forming region in Taurus, at 5.2 % us- 
ing Mgas = 1450 Mq from ^^CO (Warin et al. 1996) and 
Mstars = 79 Mq. We also estimated the SFE of the NGC 
2024 at 17 % using Mgas = 896 Mq from C^^O (Lada et 
al. 1991a,b) and the C^^O-"CO mass correlation (Ridge 
et al. 2003), and Mstars = 182 M p. Of 76 embedde d 
clusters in the solar neig hborhood (|Lada fc Ladall2003[ ). 



the SFEs for 30 embedded clusters, whose Mgas could be 
found in the literatures, are estimated in the same way, 
resulting in the SFEs widely ranging from 2.3 to 57 % 
with an median of 15 %. Median size and mass of the 
cluster-forming cores are 0.6 pc (0.32-1.0 pc in FWHM) 
and Mgas = 480 Mq (27-4000 Mq), respectively, while 
median size and mass of embedded clusters are 0.62 pc 
(0.3-1.9 pc) and Mstars = 66 Mq (25-340 Mq), respec- 
tively. 

The spatial distribution of the Cloud 2 cluster members 
(red crosses) and the molecular cloud cores (blue contour 
for ^^CO and gray scale for ^^CO) are shown in Fig. [TTl 
Mstars of the Cloud 2-N and -S clusters are estimated 
at 48 and 46 Mq, respectively, by directly counting the 
stellar mass of all cluster members larger than mass de- 
tection limit (< O.IMq) using the model KLF at the age 
of 0.5 Myr constructed in § 16.21 Even if we estimate 
Mstars with the stars down to 0.017 Mq as the same way 
in Lada & Lada (2003), the Mstars increases only ~ 5 
% . Both sizes of the clusters are typical, whose radii 
are ~ 1 pc for Cloud 2-N and ~ 0.6 pc for Cloud 2-S. 
For Mgas, because of the large distance to Cloud 2 and 
the relatively low spatial resolution of radio observations 
compared to NIR observations, the region of associated 
core must be carefully assessed. In Fig.[TT]the Cloud 2-N 
core has two sub-peaks which have similar flux, while the 
Cloud 2-S core appears to be single. Because the Cloud 
2-N cluster is associated with the northern sub-peak, we 
have estimated the Mgas for the Cloud 2-N cluster us- 
ing only the northern part of the cloud. The resultant 
Mgas for the Cloud 2-N and Cloud 2-S clusters are 360 



Mq and 310 Mq with the radii of 2.7 and 1.7 pc in 
FWHM, respectively (Saito et al. 2007). As a result, 
the SFEs of Cloud 2-N, -S clusters are estimated at 12 
and 13 %, respectively. These SFEs are comparable to 
embedded clusters in the solar neighborhood. For FOG, 
there have been s everal studi e s on S EE of star forming 
molecular clouds. ISnell et all ()2002f ) found that the ra- 
tio of far-infrared luminosity to molecular cloud mass in 
FOG is comparable to that of molecular clouds in the 
solar neighborhood, suggesting SFE in FOG is similar to 
that in the solar neighborhood. Our estimates of SFEs 
in EOG are consistent with these results for FOG. It is 
also interesting to note that the SFEs for the two inde- 
pendent clusters are similar despite the large difference 
in the stellar densities (see § 16. ip . 

Note that the estimated SFEs still have some system- 
atic uncertainties that should be carefully considered. 
Firstly, it is possible that some stars were not counted 
because of confusion with other stars. For example we 
might miss more close binaries at this large distance. 
However, even if all stars are binary with the same mass, 
the SFEs of the Cloud 2-N and -S clusters are 24 % and 
26 % respectively. This is not unusual. Secondly, it is 
possible that we are not still resolving the cloud cores 
because of the large distance to Cloud 2, thus overesti- 
mating the cloud mass by including the mass of unrelated 
cores. Actually, the radii of Cloud 2 cores (> 2 pc) are 
about twice the maximum radius of nearby cores (< 1 
pc), there is a possibility that the Cloud 2 cores arc not 
resolved by a factor of two. In this case, SFEs are un- 
derestimated. Lastly, Mgas may be underestimated given 
low metallicity of Cloud 2 because we assume the normal 
fractional abundance of X{}^CO) at 2 x 10~^ in deriv- 
ing Mgas- Since the metallicity is lower by a factor of 
5 than that in the solar neighborhood, the Xi^^CO) is 
likely 5 times larger than the assumed value. As a result, 
Mgas may become larger accordingly. In this case, SFEs 
are overestimated. Of these uncertainties the estimate 
of Mgas is most critical. For more accurate estimate of 
SFEs, radio observation with higher spatial resolution 
(e.g., with interferometry) is ncccssary and a more accu- 
rate CO-H2 conversion factor Ar(^^CO) as a function of 
metallicity is needed. 

7. CONCLUSION 

We have conducted a pilot study of star formation in 
EOG using the two embedded clusters in Cloud 2. Our 
conclusions can be summarized as follows: 

1. We have obtained deep near-infrared images of two 
embedded clusters in EOG, the Cloud 2 clusters, 
with Subaru 8.2 m telescope and the IRCS infrared 
imager. This is the first deep infrared imaging of 
star forming regions in EOG, with the mass detec- 
tion limit close to the substellar regime (< O.IA/q). 

2. Due to the faint detection-limit, we investigated 
the IMF and age of the clusters in EOG by iden- 
tifying the peak magnitude of the KLF. The im- 
derlying IMF of the cluster down to the detection 
limit is not significantly different from the typical 
IMFs in the field and in the nearby star clusters. 
For investigating the behavior of IMF (KLF) in the 
substellar regime, even deeper imaging is required 
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and we might find a different IMF than tlie solar 
neighborhood. We have estimated the age of the 
clusters to be ~ 0.5-1 Myr. Additionally, SFEs 
of the Cloud 2 clusters were comparable (~ 10 %) 
within the uncertainties to nearby embedded clus- 
ters. Despite these identical characteristics of the 
Cloud 2 clusters, the estimated stellar density sug- 
gests that different types of star formation are on- 
going: an "isolated type" star formation is ongoing 
in Cloud 2-N, while a "cluster type" star formation 
is ongoing in Cloud 2-S. 

3. We summarized the existing models of near- 
infrared characteristics of low metallicity stars to 
derive the mass-luminosity relation and color-color 
diagram for low metallicity stars. We concluded 
that dwarf track and reddening vector in JHK 
color-color diagrams do not strongly depend on the 
metallicity down to that of the stars in Cloud 2 
(~ —0.7 dex). Therefore, we did not include the 
effects of low metallicity in the following analysis 
using color-color diagram, such as the identifica- 
tion of cluster members, the estimation of extinc- 
tion and color excess. However, the effect could 
be significant for even lower metallicity and more 
theoretical studies and compilations of the near- 
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Fig. 1. — Left: JHK pseudocolor image of the Cloud 2-N cluster. North is up and east is left. The field of view is about 1' X 1'. 
The coordinate of the field center is (02000 1 feooo) = (02''48'"42.5'', +58°28'56''0) with an uncertainty of less than 0.5 arcseconds. Right: 
Spatial distribution of Cloud 2-N cluster members. Filled circles and open circles represent cluster members and field stars, respectively. 
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Fig. 2. — Left: JHK pseudocolor image of the Cloud 2-S cluster. North is to up and east is left. The field of view is about 1' X 1'. 
The coordinate of the field center is (0200O1 feooo) = (02''48™'28.7'', +58°23'34'.'6) with an uncertainty of less than 0.5 arcseconds. Right: 
Spatial distribution of the members of the Cloud 2-S cluster, the Cloud 2-S sub-cluster and fi eld stars are shown with filled c ircles, circles 
wit h cro ss, and open circles, respectively. The Cloud 2-S sub-cluster was designated at IRS3 in lKobavashi fc Tokunagal (120001 ) as discussed 
in §0 
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Fig. 3. — Le/t: (J — H) vs. {H — K) diagram in the MKO photometric system. The dwarf star track (solid line) is constructed from 
colors of dwarf stars observed with MKO filters (Hewett et al. 2006, open circle, B9V-M6; Leggett et al. 2002, plus symbol, M4V-M6V; 
Leggett et al. 2006, cross symbol, M4-M6). See the details in the main text. CTTS locus (dot-dashed line) is estimated by least-square 
fit of the CTTS data (star symbol). The reddening vector for Ay = 5 mag is shown with a dashed arrow from M6V location on the dwarf 
star track. Right: Comparison of dwarf star tracks (solid lines), CTTS loci (dot-dashed lines), and reddening vectors (dash ed arrows) i n 
MKO (black) and other photometric syste ms (gray). For the latter, the dwarf star track is in Johnson-Glass system (Bc ssel & : Bretl)[l988l) . 
the reddening vector is in Arizona system I IRieke fc Lebofskvlll985h . and the CTTS locus is in CIT system IIMever et aLlUQOTfT 




Fig. 4. — Dwarf star tracks of differing metallicit ies in JHK co lor-color diagram. Note that all the colors arc shown in CIT system. 

Color lines are low-mass star tracks from Table 6 in lLeggettI p99l '). YD ([M/H] ~ 0.0; M0-M6), OD ([M/H] 0.5; M0-M5.5), and H 

([M/H] < —1.0 ; K7-M 3) represent stars in young disk, old disk, and halo, respectively. The gray line shows the standard dwarf track from 
IBessel fc BrettI I I1988I) after converting the color system from Johnson-Glass to CIT using the transformation in lLeegettI II1992I) . 
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Fig. 5. — (J — H) vs. {H — K) color-color diagrams of the Cloud 2-N cluster (top) and the Cloud 2-S cluster (bottom). Identified cluster 
members and field stars are shown with filled circles and open circles, respectively. Cloud 2-S sub-cluster members are shown with circled 
crosses (bottom). Only stars detected with more than 5cr in all J H K-hands are plotted. The dwarf star track, CTTS locus, and reddening 
vector are based on the MKO photometric system (see Fig.|3]and discussion in §[3l §|4)l. For stars that are not detected with 5 cr in J-band, 
we showed the J — H upper limit. 




Fig. 6. — Av distribution of Cloud 2 cluster members (thick lines) and field stars (thin lines). Average Ay values of the members in 
Cloud 2-N and -S clusters are 7.2 and 6.1 mag, respectively, while those field stars in the Cloud 2-N and -S clusters arc both 2.2 mag, 
respectively. Note that Cloud 2-S sub-cluster members are not counted in the right plot. 
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Fig. 8. — Avatar VS. cluster ra dius for the Cloud 2- N cluster (open circle), Cloud 2-S cluster (filled circle), plotted along with the data for 
nearby embedded clusters from lLada fc Ladal II2003I . squares) and Carpenter (2000, stars). The solid line shows the constant stellar density 
of 50 pc~^. The gray open and filled circles show the points for the Cloud 2 clusters for the case that all the stars are unresolved binaries. 
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Fig. 9. — Top: The raw KLF (gray linos) and completeness-corrected KLF (black lines) of Cloud 2 clusters. The observed KLFs for the 
Cloud 2-N and Cloud 2-S clusters are shown with filled circle and open circle, respectively. Error bars show the uncertainties due to Poisson 
statistics. Dashed lines show the limiting magnitudes of 5 cr detection (Cloud 2-N: 20.9 mag, Cloud 2-S: 20.4 mag). Bottom: Detection 
completeness curves derived from the simulations using the artificial point sources. See the details in the main text. 
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Fig. 10. — Comparison of completeness-corrected KLF (thick line) with model KLFs of various ages (colored lines) for the Cloud 2-N 
cluster (left) and Cloud 2-S cluster (right). Red, green, blue, and orange lines represent the model KLFs of .1 Myr, 0.5 Myr, 1 Myr, and 
2 Myr, respectively. The underlying IMF for the model KLFs is the Trapezium IMF bv lMuench et al.l II2002I ). 
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Fig. 11. — Spatial distribution of Cloud 2 clusters (red erosses) and Cloud 2 eores (blue contour for ^^CO and gray scale for ^■^CO from 
Saito et al. (2007)). The Cloud 2-N has two sub-peaks, while the Cloud 2-S core appears to be single core. The Cloud 2-N cluster seems to 
be associated only with northern core of the Cloud 2-N. The contour interval is 0.05 K km s^'^ and ranges from 0.20 to 0.45 K km s'^ and 
the gray-scale range is from 0.0 to 1.7 K km s^^. The temperature scale is in T^. The beam size is about 19" for both ^-^CO and ^^CO. 



